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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

TECHNICAL  NOTE  2492 

A  METHOD  OF  SOLVING  THE  DIRECT  AND  INVERSE  PROBLEM  OF  SUPERSONIC 
FLOW  ALONG  ARBITRARY  STREAM  FILAMENTS  OP 
REVOLUTION  IN  TURBOMACHINES 
By  Chung-Hua  Wu  and  Eleanor  L.  Costilow 


SUMMARY 

I 

Analysis  of  the  supersonic  flow  in  two  two-dimensional  high- 
solidity  cascades  and  in  a  partly  supersonic  symmetrical  nozzle  shows 
that  there  is,  in  general,  significant  deviation  of  the  mean  streamline 
shape  from  that  of  the  mean  Blade  line  and  that  the  effect  of  blade 
thickness  and  blade  curvature  on  the  specific  mass  flow  along  the  mean 
streamline  is  to  increase  the  specific  mass  flow  along  the  mean  stream¬ 
line  about  9  percent  above  that  given  by  a  one-dimensional  estimate  - 
In  order  to  determine  these  effects  more  accurately  for  turbomachines 
of  arbitrary  hub  and  casing  shapes  to  be  used  for  the  three-dimensional 
through-flow  calculation,  a  method  is  developed  for  the  determination  of 
the  supersonic  flow  along  stream  surfaces  of  revolution  in  turbomachines. 
In  this  method,  the  shapes  of  the  stream  surfaces  are  arbitrary,  and  the 
method  also  takes  into  account  the  distance  between  adjacent  stream 
surfaces,  which  varies  along  the  flow  path.  Thus,  the  method  can  be 
applied  to  turbomachines  with  arbitrary  hub  and  casing  shapes. 

In  addition  to  their  use  for  direct  problems,  these  equations  can 
be  used  to  design  blade  elements  in  supersonic  flow  along  an  arbitrary 
stream  filament  of  revolution  in  turbomachines . 


INTRODUCTION 

Recent  investigations  of  the  applicability  of  supersonic  flow  in 
compressors  have  shown  the  desirability  of  using  such  flow  to  increase 
the  pressure  ratio  per  stage  (references  1  and  2) .  For  axial-flow 
compressors  having  blades  with  short  radial  length,  the  over-all 
performance  analysis  is  often  based  on  a  one- dimensional  approximation, 
in  which  only  average  values  in  the  channel  are  considered  (for  example, 
references  3  and  4) .  For  radially  long  blades  and  variable  root  or  tip 
radii,  methods  are  proposed  for  analyzing  the  flow  by  a  three- 
dimensional  "through-flow"  calculation,  which  considers  the  axial  and 
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radial  variations  of  the  flow  hut  only  a  mean  value  in  the  circumferen-  « 

tial  direction  (references  5  and  6) .  The  flow  on  the  relative  mean 
stream  surface,  which  divides  circumferentially  the  mass  flow  passing 
through  the  channel  between  two  blades  into  two  equal  parts,  is  taken 
to  represent  the  mean  flow  through  the  blading.  It  is  suggested  in 
reference  5  that  the  shape  of  this  mean  stream  surface  and  the  correction 
factor  b  for  a  finite  number  of  thick  blades  be  obtained  from  the 
analysis  of  a  number  of  two-dimensional  flows  along  general  surfaces  of 
revolution  starting  at  different  inlet  radii.  A  method  was  therefore 
developed  at  the  MCA  Lewis  laboratory  to  determine  the  flow  variation  on 
these  flow  surfaces  in  a  supersonic  turbomachine.  When  the  flow  surfaces 
may  be  approximated  by  cylindrical  surfaces,  the  flow  equations  reduce  to 
the  usual  plane  flow  where  the  hodograph  characteristics  are  applicable 
for  the  flow  analysis  if  the  entrance  shock  is  weak. 

In  order  to  illustrate  the  effects  of  blade  thickness  and  blade 
curvature  and  to  determine,  in  general,  the  correlation  between  the 
shapes  of  the  mean  streamline  and  the  mean  blade  line,  flow  on  the  mean 
streamlines  was  determined  for  two  supersonic  cascades  and  a  partly 
supersonic  hyperbolic  nozzle. 

This  method,  in  addition  to  its  use  for  flow  analysis  of  a  given 
blading,  can  be  used  to  design  blade  elements  along  stream  filaments  of 
revolution  having  arbitrary  thickness  variation  along  the  flow  path. 


SYMBOLS 

The  following  symbols  were  used  in  this  report: 

A  area  normal  to  velocity 

a  local  speed  of  sound 

b  correction  factor  for  finite  number  of  thick  blades 

c  line  tangent  to  characteristic  curve 

Cp  specific  heat  at  constant  pressure 

h  static  enthalpy 

I  h  +  |w2  -  |o)2r2 

J,K,L  coefficients  of  partial  derivatives  of  \|r 
M,N 
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orthogonal  coordinates  on  mean  surface  of  revolution 


local  Mach  number 


mass  flow  “between  suction  surface  of  “blade  and  streamline, 


pw  a| 


“blade  pitch  or  spacing 
gas  constant 


radial  distance  from  axis  of  machine  (fig.  l) 
entropy 


k 


static  temperature 

“blade  thickness  in  the  y-direction 

resultant  relative  velocity 

rectangular  coordinates  for  cylindrical  flow  with  y-axis* 
chosen  along  line  joining  leading  edges  of  blade 

W 

-1  CD 

relative  flow  angle,  tan 

1 1 

ratio  of  specific  heats 

slope  of  characteristics  on  mean  surface  of  revolution 
Mach  angle,  sin--3-  ^ 

distance  measured  from  suction  to  pressure  surface  normal  to 
resultant  velocity 

mass  density 

slope  of  mean  surface  of  revolution  in  meridional  plane, 

W 

tan--3-  rr~ 


normal  thickness  of  stream  filament  of  revolution 


stream  function 
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a)  angular  velocity  of  blade 

Subscripts: 

i  inlet  condition 

2,cp  meridional  and  circumferential  components 

m  on  mean  streamline 

r,z  components  in  r-  and  z-directions,  respectively 

s  initial  supersonic  state 

t  at  throat  of  nozzle 

x,y  components  along  x-  and  y-directions,  respectively 

1,2  first  and  second  family  of  characteristics,  respectively 


METHOD 

On  a  general  surface  of  revolution  defined  by  the  orthogonal 
coordinates  l  and  cp  (figs.  l(a)  to  1(c)),  the  supersonic  flow  of  a 
fluid  along  a  stream  filament  of  revolution  of  varying  normal  thickness 
T  =  t(Z)  is  described  by  the  following  forms  of  the  flow  equations .  By 
analyzing  the  flow  going  in  and  out  of  the  element  shown  in  figure  l(d), 
the  continuity  equation  for  steady  flow  may  be  expressed  by 


S(TpW?r)  d(TPw  ) 

ST-  + 

When  equation  (l)  is  expanded  and  the  relation 


(1) 


d  In  p 


is  used,  equation  (l)  becomes 


ST  <WZTr) 


(TWj  +  -L 


■cp 


a1 


Sh 

ST 


S 


(2) 
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For  adiabatic  flow,  the  entropy  along  a  streamline  remains  constant  for 
a  nonviscous  fluid;  that  is 


Ds  TI  5s  ,  wcp  5s 
Dt  3T  +  r  5V  0 

which  allows  the  last  term  to  he  dropped  from  equation  (2) 
perfect  gas, 


For  a 


h  =  cpT 


a2  =  (r  -  l)  h 

When  this  substitution  is  made  for  a2  and  the  expression  is  multiplied 
1 

by  iF-1,  the  equation  of  continuity  (2)  becomes: 

^(TW2rhT"1)+  -  0  (3) 

Expressed  in  the  coordinates  l  and  qp,  the  equation  of  steady  motion  of 
a  nonviscous  fluid  as  given  by  equation  (l4a)  in  reference  5  in  the 
circumferential  direction,  is 

I  If  '  Ir2  -(t2+  2a))sin  a  -  ^  (isi-  0  (4) 

From  equation  (3) ,  a  stream  function  is  defined  by: 

1 

H  =  ^TWjh7'1  (5a) 


||  =  -  T^hr  1  (5b) 

Substituting  equations  (5)  into  equation  (4)  results  in  the  follow¬ 
ing  relation: 
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i  ft 

V  a 2  /  52^  a2  r  dcp  dl  \  a2  /  r2 


1_ 

2 


sin 


cr  + 


a2  St 

TST  + 


Si 

ST 


S\|/ 

ST' 


A 


i|i  =  0 

r  dcp 

(6) 

If  the  symbols  J,  K,  L,  M,  and  N  are  used  to  represent  the  coefficients 
of  the  partial  derivatives  in  equation  (6),  it  may  he  rewritten  as 


2ao 


sin  a 
W, 


T  5s 

r  Sep 


9 

a^r 


5l 

5cp 


52^ 

522 


+ 


2K  52\|r  L  52\|r  ,  M  5^  N  5^  _ 

r  5cp  dl  r2  ^p2  5Z  r  £cp 


0 


The  characteristics  of  equation  (7)  are  then 


(?) 


J 


2K 


(8) 


For  equation  (8),  two  real  characteristics  exist  for  supersonic 
flow,  the  slopes  along  which  are 


When  the  polar  coordinates 

^  =  sln"1(i)  =  sin_1(t)  =  sin_1 


acp\  _  k  Vk2-jl 
a  i]1  j  "  j 


acp\  k  .  7k2-jl 
ai/2  "  J  +  J 

transformation  that 


-  is  used, 

A 2  +  V 


(9a) 

(9b) 

,  Wcp 

p  =  tan-J-  —  ana 
WZ 

the  relation  of 


to  the  flow  angle  ana  the  Mach  angle  is 
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=  tan  (p  +  |i) 


(lOa) 


A2  =  tan  (p  -  (i) 


(lOt) 


By  the  use  of  equations  (7)  and  (9),  the  rate  of  change  of 
^-derivatives  along  the  characteristics  with  respect  to  Z  is  obtained 
(see  reference  5  for  procedure) :  Along  c^, 


and  along  c2 


d  ( S\|r  \  ^2  d_  ^  M  Si|r  _N_ 

dZ  V  ST/  r  dZ  VScp/  J  SI  +  Jr  Sep  =  0 


d/’Si|r'\  ^ld/"S'^)iMS\lr  Nch|r 

iTv^ry  ~dTvs^/  +  jsr  +j?s^=0 


(lla) 


(in) 


For  ease  in  computation,,  it  is  more  convenient  to  express  equa¬ 
tions  (id  in  terms  of  the  magnitude  of  the  resultant  velocity  W  and 
the  flow  angle  P  as  follows:  Substituting  equations  (5)  in  equa¬ 
tions  (ll)  and  differentiating  yields  the  following  expression  along  Cj_: 


/-v  x  t  dh  /  dw7  dW„\ 

(*2  wz  -  V  fs  it  +t^2  dr  -  dry+  (A 


2  wz  -  V  ir + 


1  (NWZ 


-  MWjp)  +  TWZ  sin  cr  =  0 


Similarly,  along  c2, 


w7  -  vj  II  + 

1  dZ 


J  (HWj  -  Wq)  + 


TW 7  sin  a=0 
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With. 


W^  =  W  cos  3 
=  W  sin  3 


and 


ah 

d  l 


W2  a>2r2  \ 
2  2  J 


dl 
d  l 


it  aW  2 

W  *rr—  +  a)  r  sm  a 
d  7 


and  with  the  use  of  equations  (10),  equations  (ll)  become 


1  dW 
W  dl 


dfi  ? 

-  tan  \i  g  -  tanc  p, 


cd  r  sin  a 


1_  dl 
n2  dl 


Xg  sin  a  cos  p 


1  £[ 

+  T  dl  +  i^Xg  cos  0  -  sin  3) 


1  N  cos  p  -  M  sin  P 
J  (X2  cos  p  -  sin  p) 


=  0 


(12a) 

i  —  +  tan  (i  ^2.  -  tan2  p 
W  dl  *  dl  K 

(12b) 


orr  sin  a  +  1  dl 
Q2  n2  dl 


+  15!  + 

T  dl 


Xn  sin  a  cos 


+  1  N  cos  p  -  M  sin  p 


r(X-^  cos  p  -  sin  p)  J  (X^  cos  p  -  sin  p) 


=  0 


along  C]_  and  c 2,  respectively. 

From  two  points  a  and  b  (fig.  2)  where  the  flow  is  known,  the 
tangents  to  the  characteristics  curves  at  these  points  are  determined 
from  equations  (9a)  and  (9b),  respectively,  and  intersect  at  a  point  c. 
By  writing  equation  (l2a)  (equation  along  cj)  at  a  and  equation  (l2b) 
(equation  along  eg)  at  b  in  finite  difference  form  and  solving  simul¬ 
taneously,  the  new  flow  values  at  c  are  determined;  thus,  along  c^ 
at  a, 


1 

Wa 


tan  na 


co^r  sm  a 


tan^  [ia 


and  along  c2  at  b, 


(13) 


tan 


Pc  -  Pb 

lc  '  2h 


tan2 


9 

ay  r  sm  o-^ 

^b  - n -  + 

er 


(14) 


The  unknowns  in  both  equations  (13)  and  (l4)  are  Wc  and  |3  ;  A 7  is 
determined  by  the  intersection  of  the  tangents  to  the  characteristics 
from  points  a  and  b  so  that  7C  -  7a  and  7C  -  7 ^  are  known  in  the 
equation. 
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For  a  characteristic  curve  that  intersects  the  passage  boundary, 
such  as  is  shown  at  d  or  e  in  figure  2,  the  situation  is  a  little 
different.  In  a  direct  problem  with  a  given  blade  configuration,  the 
angle  (3  at  the  intersecting  boundary  point  is  known.  In  the  inverse 
or  blade-design  problem,  either  (3  at  the  intersecting  boundary  point 
is  known  from  a  desirable  turning  specified  at  that  point  or  the  magni¬ 
tude  of  the  resultant  velocity  W  at  the  intersecting  boundary  point  is 
known  from  a  specified  velocity  distribution  on  the  blade.  In  all  cases, 
either  equation  (13)  or  (14)  written  for  the  characteristic  from 
d  or  e  determines  the  unknown  W  or  (3  at  the  intersecting  boundary 
point,  whereas,  as  shown  previously,  two  simultaneous  equations  along 
intersecting  characteristics  are  required  inside  the  passage.  In 
general,  after  W  and  (3  are  determined  at  a  new  point,  other  fluid 
properties  at  that  point  are  calculated  and  the  characteristic  directions 
constructed  to  repeat  the  process  progressively  downstream.  (For  irrota- 
tional  inlet  flow  with  uniform  I  and  s  and  when  the  rotationality 
introduced  by  the  entrance  shocks  is  neglected,  I  and  s  are  uniform 
throughout  for  adiabatic  flow,  and  the  derivatives  involved  in  equa¬ 
tions  (12)  become  zero.  Otherwise,  the  derivatives  are  to  be  evaluated 
by  using  the  values  at  the  inlet  on  different  streamlines,  the  constancy 
along  the  streamline  before  and  after  shock,  and  the  changes  across  the 
shock.) 

Thus  a  method  is  available  either  to  analyze  the  blade-to-blade 
flow  variation  along  a  given  arbitrary  stream  filament  of  revolution  in 
a  supersonic  turbomachine,  or  to  design  the  blade  element  on  a  given 
arbitrary  stream  filament  of  revolution  for  a  specified  turning  distri¬ 
bution  or  a  specified  velocity  distribution  around  the  blade.  The  con¬ 
figuration  of  the  stream  filaments  (fig.  1(a))  is  taken  from  a 
through-flow  calculation  (references  5  and  6) .  For  a  direct  problem, 
successive  calculations  between  the  through-flow  calculation  and  the 
present  calculation  are  necessary  until  the  solution  converges.  The 
through- flow  calculation  gives  the  configuration  of  stream  filaments  of 
revolution,  and  the  calculation  on  these  stream  filaments  of  revolution 
gives  the  shape  of  the  mean  stream  surface  and  the  factor  b  used  in 
reference  5,  which  accounts  for  the  effect  of  blade  thickness  and 
curvature.  In  am  inverse  problem,  the  calculation  is  shortened  if  the 
estimated  values  of  the  factor  b  used  in  the  through-flow  calculation 
for  a  desirable  blade-thickness  distribution  results  in  a  blade  dimension 
which  has  a  thickness  distribution  close  enough  to  the  desirable  one. 

In  order  to  give  a  general  idea  of  the  order  of  magnitude  of  this 
factor  b  for  the  effect  of  the  blade  thickness  and  curvature  on  the 
mean  flow  and  the  closeness  of  the  shape  of  the  mean  stream  surface  to 
the  blade  mean  surface,  an  analysis  is  made  for  the  simpler  case  of 
cylindrical  flow  utilizing  available  data  for  two  two-dimensional 
high- solidity  cascade  recently  investigated  by  Liccini  at  the  NACA 
Langley  laboratory  (unpublished  data)  and  for  a  partly  supersonic  nozzle 
investigated  by  Emmons  (reference  7) . 
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ANALYSIS  OF  FLOW  IN  TWO  SUPERSONIC  CASCADES  AND 

ONE  PARTLY  SUPERSONIC  NOZZLE  • 

The  shapes  of  the  two  high- solidity  90°  turning- angle  cascades  are 
shown  in  figure  3.  Each  cascade  has  an  inlet  Mach  number  of  1.78  and 
a  10°  wedge  angle  at  the  leading  edge.  (The  blade  is  so  designed  that 
the  leading-edge  shock  is  cancelled  after  the  first  reflection.)  The 
shape  of  the  partly  supersonic  nozzle  is  shown  in  figure  4.  (Because 
the  nozzle  is  symmetrical,  only  half  of  it  is  shown.) 

In  the  case  of  the  symmetrical  nozzle,  the  center  line  is  the  mean 
streamline.  For  the  two  90°  cascades,  the  mean  streamline  (the  one 
which  divides  the  mass  flow  in  the  channel  into  two  equal  parts)  has  to 
he  established  from  the  data  given  on  the  characteristic  net  of  the 
blade  design.  The  mass  flow  m  at  any  point  in  the  channel  is  deter¬ 
mined  by  integrating  across  the  channel  the  local  specific  mass  flow 

rt 

given  by  the  characteristic  net;  thus,  m  =  J  pW  d|.  The  values  of 

°  o 

the  local  specific  mass  flow  can  he  obtained  by  using  the  table  given  in 
reference  8  or  similar  tables.  The  shape  of  the  mean  streamline  is  to 
be  compared  with  that  of  the  mean  blade  line,  which  is  obtained  by 
taking  the  mean  of  the  blade  coordinates  in  the  pitch  direction. 

The  value  of  the  specific  mass  flow  W  or  Wx  is  then  obtained 
along  the  mean  streamline  downstream  of  its  intersection  with  the 
leading-edge  shock.  The  discontinuity  in  the  flow  curves  at  the  point 
where  the  mean  streamline  intersects  the  "reflection"  is  smoothed  out 
in  order  to  emphasize  the  general  trend.  These  values  are  to  be  com¬ 
pared  with  the  corresponding  variation  in  the  channel  width  (or  area) 
perpendicular  to  W  or  Wx  in  order  to  determine  if  there  is  any  corre¬ 
lation  between  the  variations  in  specific  mass  flow  on  the  mean  stream¬ 
line  and  the  corresponding  area. 

In  the  case  of  the  nozzle,  the  flow  variation  on  the  mean  stream¬ 
line  is  obtained  from  the  Mach  number  given  in  reference  7.  The  vari¬ 
ation  in  specific  mass  flow  is  compared  with  the  channel- area  ratio. 

The  results  obtained  in  the  two  90°  turning  passages  are  given  in 
figures  3  and  5  to  9 .  The  results  obtained  with  the  symmetrical  nozzle 
are  shown  in  figures  4,  10,  and  11. 

In  figures  3(a)  and  3(b),  the  mean  streamline  is  compared  with  the 
mean  channel  line.  In  figure  3(a),  where  the  cascade  is  of  the  higher 
solidity  and  has  thinner  blades,  the  shape  of  the  mean  streamline 
closely  approximated  that  of  the  mean  camber  line;  however,  in  the 
second  cascade  (fig.  3(b)),  which  has  a  lower  solidity  and  thicker 
blades,  the  mean  streamline  deviates  more  from  the  mean  blade  line. 
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The  specific  mass  flow  on  the  mean  streamline  based  on  the  resultant 
velocity  as  shown  in  figure  5  indicates  clearly  that  the  effect  of  blade 
thickness  and  curvature  on  the  through  flow  is  to  increase,  in  both 
cases  investigated,  the  specific  mass  flow  along  the  mean  streamline  a 
chordwise  average  of  9  percent  over  that  given  by  the  area  reduction 
(due  to  blade  thickness)  based  on  a  one- dimensional  calculation.  If  a 
comparison  of  specific  mass  flow  in  the  axial  direction  with  variation 
of  channel  width  in  the  y-direction  is  made  (fig.  6),  as  is  done  in 
reference  9  for  subsonic  flow,  corresponding  trends  are  obtained,  and 
the  increase  in  specific  mass  flow  again  averages  about  9  percent 
chordwise  over  the  area  ratio.  This  increase  is  more  than  two  times  the 
increase  obtained  in  reference  9. 

The  variations  of  axial  and  tangential  velocity  components  on  the 
mean  streamline  are  shown  in  figures  7  and  8,  respectively.  The  vari¬ 
ations  of  the  velocity  across  the  channel  in  the  two  cascades  at  the 
stations  indicated  on  figure  3  are  shown  in  figure  9. 

The  result  obtained  in  the  nozzle  problem  is  shown  in  figures  10 
and  11.  The  results  shown  in  figure  10  correspond  to  the  nozzle  shown 
in  figure  15  of  reference  7,  which  includes  a  small  portion  of  super¬ 
sonic  flow  in  the  divergent  portion  of  the  nozzle;  whereas  figure  11 
corresponds  to  figure  16  of  reference  7,  with  supersonic  flow  extending 
throughout  the  divergent  portion  and  the  Mach  number  increasing 
beyond  2.  For  the  subsonic  flow  in  the  convergent  portion,  the  local 
specific  mass  flow  and  area  are  normalized  by  the  values  at  the  throat; 
for  the  flow  in  the  divergent  portion  (partly  supersonic,  fig.  8,  and 
purely  supersonic,  fig.  9),  the  local  specific  mass  flow  and  area  are 
normalized  with  respect  to  the  conditions  where  the  flow  is  initially 
supersonic.  The  comparison  between  the  specific  mass  flow  and  area 
ratio  is  similar  to  the  previous  result.  The  increase  in  the  specific 
mass  flow  on  the  mean  streamline  is  entirely  due  to  the  area  reduction 
(no  turning)  and  is  about  8  percent  higher  than  a  one- dimensional 
correction. 

The  results  obtained  in  these  analyses  give  some  idea  regarding  the 
closeness  between  the  shape  of  the  blade  mean  line  and  mean  streamline 
and  of  the  order  of  magnitude  of  the  blade  thickness  and  curvature 
effect  for  through-flow  calculations.  For  an  approximate  through- flow 
solution  of  a  similar  blading,  these  results  may  be  directly  applied. 

For  different  types  of  blading,  similar  analyses  can  be  made  for  the 
two-dimensional  flow  on  cylindrical  surfaces  or  surfaces  of  revolution, 
which  are  formed  by  fluid  particles  initially  lying  on  circular  arcs  at 
a  number  of  radii  upstream  of  the  machine  by  using  the  general  relations 
derived  herein.  The  mean  streamline  obtained  at  the  different  radii  can 
be  joined  together  to  form  the  mean  stream  surface.  The  ratio  of 
(p^x^m  "to  (pW^i  obtained  in  these  calculations  can  be  taken  as  the 
thickness  correction  factor  b  in  the  relations  between  stream  function 
and  velocities  (equation  (ill)  of  reference  5). 
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More  knowledge  of  this  kind  for  a  number  of  typical  bladings  will 
be  very  useful  in  the  design  of  these  bladings.  From  this  knowledge  and 
a  certain  blade-thickness  distribution,  which  is  desired  from  the  blade 
strength  and  other  considerations,  a  good  estimate  of  the  factor  b  can 
be  made  and  used  in  the  through-flow  calculation.  From  this  solution, 
the  blade  section  on  a  number  of  surfaces  of  revolution  can  be  designed 
by  the  method  given  in  the  report. 


SUMMARY  OF  RESULTS 

Analysis  of  the  supersonic  flow  in  two  two-dimensional  90°  turning 
passages  of  high  solidity  and  in  a  partly  supersonic  symmetrical  nozzle 
shows  that  there  is,  in  general,  significant  deviation  of  the  mean 
streamline  shape  from  that  of  the  mean  blade  line  and  that  the  effect  of 
blade  thickness  and  blade  curvature  on  the  specific  mass  flow  along  the 
mean  streamline  is  to  increase  the  specific  mass  flow  along  the  mean 
streamline  about  9  percent  above  that  given  by  a  one- dimensional 
estimate.  In  order  to  determine  these  effects  more  accurately  for 
turbomachines  of  arbitrary  hub  and  casing  shapes  to  be  used  for  the 
three-dimensional  through- flow  calculations,  a  method  is  developed  for 
the  determination  of  the  supersonic  flow  along  stream  surfaces  of  revo¬ 
lution  in  turbomachines.  In  this  method  the  shapes  of  the  stream 
surfaces  are  arbitrary,  and  the  method  also  takes  into  account  the 
distance  between  adjacent  stream  surfaces,  which  distance  varies  along 
the  flow  path.  Thus,  the  method  can  be  applied  to  turbomachines  with 
arbitrary  hub  and  casing  shapes. 

In  addition  to  their  use  for  direct  problems,  these  equations  can 
be  used  to  design  blade  elements  in  supersonic  flow  along  an  arbitrary 
stream  filament  of  revolution  in  turbomachines . 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  July  13,  1951 
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(a)  Case  a.  (b)  Case  b. 

Figure  5.  -  Comparison  "between  variations  of  specific  mass  flow  on  mean  streamline 

and  normal  channel  area  ratio. 
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Figure  8.  -  Variation  of  tangential  velocity  on 
mean  streamline . 
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Figure  10.  -  Comparison  of  variation  in  mean  specific  mass  flow  and  area  ratio 
with  short  portion  of  supersonic  flow  in  divergent  section. 


Ratio 


Figure  11.  -  Comparison  of  variation  in  mean  specific  mass  flow  and  area  ratio 
of  nozzle  with  supersonic  flow  throughout  divergent  section. 
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